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Abstract: The time-average method currently available is limited to analyzing the specific 
performance of the automatic gain control-proportional and integral (AGC-PI) based 
velocity-controlled closed-loop in a micro-electro-mechanical systems (MEMS) vibratory 
gyroscope, since it is hard to solve nonlinear functions in the time domain when the control 
loop reaches to 3rd order. In this paper, we propose a linearization design approach to 
overcome this limitation by establishing a 3rd order linear model of the control loop 
and transferring the analysis to the frequency domain. Order reduction is applied on the 
built linear model's transfer function by constructing a zero-pole doublet, and therefore 
mathematical expression of each control loop's performance specification is obtained. 
Then an optimization methodology is summarized, which reveals that a robust, stable and 
swift control loop can be achieved by carefully selecting the system parameters following a 
priority order. Closed-loop drive circuits are designed and implemented using 0.35 |iim 
complementary metal oxide semiconductor (CMOS) process, and experiments carried out 
on a gyroscope prototype verify the optimization methodology that an optimized stability 
of the control loop can be achieved by constructing the zero-pole doublet, and disturbance 
rejection capability (D.R.C) of the control loop can be improved by increasing the 
integral term. 
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1. Introduction 

Due to their wide applications in inertial navigation, automotive stability control and robots, MEMS 
gyroscopes have drawn tremendous attention of researchers in both academia and industry [1-4]. Most 
MEMS gyroscopes are based on the Coriolis force effect. The Coriolis force is proportional to the 
angular rate only on the premise of a constant vibration velocity in the drive axis. As a result, 
automatic gain control (AGC) based velocity-controlled closed-loops are commonly built to keep 
the primary resonator vibrating at its resonant frequency with constant amplitude [5-12]. The 
velocity-controlled closed-loop always introduces stability problems, which are the primary concern in 
its design process. 

The time-average method is widely applied in the stability analysis of control loops [7,9,13-15]. In 
the time-average method, the controlled amplitude is considered as a constant in one vibratory period 
to simplify the nonlinear equations of the control loop, since the variation rates of the controlled 
amplitude are usually much slower than the vibratory rate. In [9,13], the nonlinear equations are 
simplified by using the time-average method and the effects of AGC parameters on the response of the 
linearized dynamics are revealed. Similar work is also done in [14,15], which use a Lyapunov function 
to obtain a good transient response. However, in all these works, no quantitative stability criterion is 
given. In [7], a nonlinear mathematical model in time domain of a AGC-PI closed-loop is linearized by 
firstly using a time-average method. Then the mathematical expression of the stability criterion is 
derived by applying the Routh-Hurwitz criterion to the characteristic equations. Specific transient 
response performance is discussed when the control loop is a 2nd order system with only the 
proportional controller. In reality, in order to improve the loop control accuracy, an integral function is 
imported into the controller of the closed-loop, and the AGC-PI closed-loop structure is widely used in 
the driving circuit of MEMS vibratory gyroscopes which is a 3rd order system. Nevertheless, the 
specific transient response performance of a 3rd order system can hardly be attained using the 
time-average method due to the difficulties in solving the complicated nonlinear functions [7]. 
In summary, the time-average method has two main drawbacks: firstly, the entire stability analysis 
process is tedious and exhausting as it solves nonlinear functions "in the time domain"; secondly, it 
can not be used to analysis the specific performance of the 3rd order controlled loop. 

To overcome these two drawbacks of the time-average method mentioned above, a linearization 
design approach is presented in this paper. A fully linear system model of the 3rd order closed-loop 
is established and thereafter the stability and performance are analyzed in the frequency domain. 
Compared to the conventional time-average method, the proposed linearization design approach has 
the following advantages: first, the stability criterion of the closed- loop can be simply obtained 
by applying a zero-pole method on the linear model; second, the mathematical expression of each 
performance specification is obtained by applying order reduction on the 3rd order linear model, so 
that the dilemma in [7] mentioned above is solved; third, an optimization methodology is summarized, 
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which reveals that a robust, stable and swift control loop can be achieved by carefully selecting the 
system parameters following a priority order. The proposed optimization methodology is verified by 
both numerical simulations and experiments. 

This paper is organized as follows: Section 2 presents the overall linearization design approach, 
including building the linear model, revealing the stability criterion, analyzing the control loop's 
performance through mathematical expressions and numerical simulation results, and summarizing an 
optimization methodology. In Section 3, the implemented AGC-PI based closed-loop drive circuits are 
introduced briefly and detailed experimental results on a gyroscope prototype are presented to verify 
the proposed optimization methodology. Finally, a conclusion is given in Section 4. 

2. Linearization Design Approach 

The topology of the AGC-PI based velocity-controlled closed-loop for the drive mode of MEMS 
gyroscope is shown in Figure 1. The displacement of sense comb fingers induces an alternating 
current, which is converted to voltage by the trans-impedance amplifier (TIA). The voltage is fed 
into the variable gain amplifier (VGA) to generate the excitation signal. The TIA, VGA and 
primary resonator of the gyroscope form an electromechanical oscillator. The vibratory velocity of this 
electromechanical oscillator is controlled to a target value V re fby an AGC, which consists of rectifier, 
low pass filter (LPF), proportional and integral (PI) controller and VGA. In the AGC, the rectifier and 
LPF extract the amplitude information of the vibratory velocity and then compare it with the target 
value V re f to obtain an error signal V error- The error signal is amplified and integrated by the PI 
controller and then transmitted to adjust the gain of the VGA. The control loop achieves a steady state 
until the error signal is zero. 

Figure 1. AGC-PI based velocity-controlled closed-loop for the drive mode of MEMS 
vibratory gyroscope. 



Resonator 




Rectifier 



LPF +(+) 



'ref 



AGC 



PI controller 



The kinetic equation of the drive mode in gyroscope is described as: 



• ^ • 2 
X+ — X+GTX- 



Q 



(i) 



where x is the displacement of the resonator, co x is the natural frequency, Q is the quality factor, m x is 
the mass and F ext is the driving force. 



Sensors 2013, 13 



12567 



2.1. Linear System Model 

The velocity-controlled closed-loop is a nonlinear system, which processes the amplitude information 
of the sinusoidal resonance signal. It is hard to analyze the behavior of a nonlinear system, especially 
when the system usually reaches a 3rd order. Fortunately, a 3rd order linear model of the closed-loop 
can be built by linearizing two nonlinear modules in the loop, as shown in Figure 2. In the figure, K C /x 
and K F /v model the gain from displacement to capacitance and the gain from the squared excitation 
voltage to the force, respectively. Ky/c models the gain of capacitance-to-voltage in TIA and K VGA 
represents the gain coefficient of VGA. In the amplitude information control path, a rectifier detects 
the envelope of the resonance signal with gain K rech and then the LPF and PI controller process the 
amplitude signal with their correspondence transfer functions. 



Figure 2. 3rd order linear system model of the velocity-controlled closed-loop. 
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The nonlinearity of the closed-loop generates arises from two modules: the primary resonator and 
the VGA. As the phase perturbation term is negligible compared with co x /(2Q) in most MEMS 
vibratory gyroscopes, Equation (1) can be reduced to a first-order system [5], expressed as: 



1 



ext _amp 



2k s + coJIQ 



(2) 



where F ext amp and X amp is the amplitude of F ext and x respectively, k is the spring constant. The 
nonlinearity of the VGA comes from the multiplication of its input and the gain control signal, both of 
which reflect the amplitude information of the vibration velocity. When the control loop gain is high 
and the perturbation term is small compared to the velocity amplitude, the error signal is close to zero 
and the magnitude of the VGA input V in _ amp can be approximated as V re /K rect [8]. As a result, with 
respect to the amplitude information control, the electromechanical resonating loop is broken at the 
VGA's input and only a linear velocity-controlled loop exists. 
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The control loop is a negative feedback system. As shown in Figure 2, defining the small signal ac 
component of the reference voltage V re f as input V in and the amplitude of the vibration velocity as 
output Vouu the closed-loop transfer function can be illustrated as: 



out 



H loop {s) 



V in l + H !oop (s) 



(3) 



where the feedback factor is 1, and the open-loop transfer function of the linearized control loop can 
be expressed as: 



H loop( S ) = K p /V 2 K C/X K V/C K VGA V refK 



CO, 



1 



G), 



Ipf 



PIV 1 CIX ± ^VIC ± ^VGA y ref r dc 



2k s + co x l(2Q) s + co t 



S + T 



(4) 



where Vd c is the DC voltage across the drive comb fingers, coi p f is the pole of LPF, K P is the 
proportional term and Kj is the integral term in the controller, x is the loss of integrator. 

The built linear system model shown in Figure 2 can be extended to arbitrary higher-order control 
loops by obeying the following rule: calculating all the significant zeros and poles into the linear 
model, which usually are located within ten times the control loop's bandwidth. 



2.2. Stability Criterion 

Figure 3 shows the three poles and the single negative zero in Equation (4). Based on the zero-pole 
method, as the phase margin compensation characteristic of negative zero, the zero should be located 
within the largest pole among the three to promise a stable loop. Therefore, the stability criterion of the 
control loop can be simply written as: 



K I / X 

— < max(^ z nf , — , t) 
K p lpf 2Q 



(5) 



Figure 3. Zero-pole map of the 3rd order velocity-controlled closed-loop. 



0.5 



^ 0 
c 

CO 

-0.5 



■I 



Pole-Zero Map 


1 1 


1 1 1 1 


1 _ j r " 




~ " L 

% f ~T~~—- 

100 \/ 80 60 




f 


K, 1 K p 


t " 




1 1 





00 -600 -500 -400 -300 -200 -100 
Real Axis 



Sensors 2013, 13 



12569 



2.3. Optimization Methodology 

Specifications including loop gain, bandwidth and phase margin are commonly used to evaluate the 
performances of a negative feedback system in the frequency domain. Mapping these specifications to 
the provided control loop, their influence on the control loop's performances is illustrated as follows. 
The loop gain determines the disturbance rejection capability of the control loop against the parameter 
variance in both mechanical structures and circuits. The bandwidth determines the recovery speed of 
the control loop and the phase margin reflects the stability of the control loop. These specifications 
correspond with transient response specifications in the time domain in [7], referred to as setting time, 
rising time and overshooting. The mathematical expressions of these three performance specifications 
are derived as follows: firstly, the loop gain of the control loop can be obtained by substituting s = 0 in 
Equation (4), expressed as: 

( f A 

(6) 



H loop( 0 )- K total-^- 



V T J 



where: 

K tota i — K F/v2 K c/x K v/c K VGA V ref V dc - -^j- (7) 

The expressions of the bandwidth and phase margin can be calculated by solving a third-order 
equation generated from Equation (4). However, the solution process is very complex. Let the negative 
zero Kj/Kp equal to the pole co x /(2Q). As a result, a zero-pole doublet is constructed and the 3rd order 
system is reduced to a second-order system. The validity of the order reduction will be revealed in 
following numerical simulations. The order-reduced transfer function from Equation (4) can be 
rewritten as: 

H loop {s) = K total -^.^- (8) 

S + Q) lpf S v ' 

in which the integrator is considered as ideal and its loss is ignored. Then the expression of bandwidth 
a phase margin can be calculated as: 

BW _ %fU X + AK ^ K P l0) w- X ) (9) 



and: 

\jl + 4K U!iul K p /m !>/ -1 

(10) 
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According to Equation (6), Equations (9) and (10), each performance specification variation 
trend versus the primary system parameters are summarized in Table 1. As shown in the table, the 
performance specifications of the control loop contradict each other. Increasing the loop gain improves 
the disturbance rejection capability and the control accuracy of the system, while sacrificing the phase 
margin as well as the stability. Extending the bandwidth increases the response speed of the control 
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system, but may reduce the phase margin and cause an unstable oscillation. As a result, some trade- 
offs should be made to achieve a precise, stable and swift control loop. 

Table 1. Performance specifications of the control loop versus primary system parameters. 



Parameters 


Loop Gain 


Bandwidth 


Phase Margin 




T 


t 


i 


K, t 


t 








— > 


t 


t 




t 






K mal T 


T 


t 


i 



T Increase; nI Decrease; — > no change; —unknown. 



The order of the consideration on the performance specification optimization is usually 
loop gain > phase margin > bandwidth. The reason can be illustrated as follows: because the 
mechanical parameter variations are dominant in the control loop and most parameters of the primary 
resonator, like quality factor Q, fluctuate slowly but with a wide range [16]. To suppress these 
mechanical parameter variations and maintain a constant vibration, the loop gain of the control loop 
should be considered first. When the gyroscope is under an external disturbance, a low phase margin 
of the control loop would introduce a large overshoot and cause a large fluctuation at the sense part of 
the gyroscope. As a result, the phase margin should be considered secondly and the bandwidth is 
the last. 

Numerical simulation results of both the 3rd order control loop linear model and corresponding 
nonlinear model are provided in the following paragraphs for three reasons: first, it is intended to 
reveal the validity of the order-reduction; second, it is used to verify the derived results in Table 1; 
third, it is used to supply the effects of parameters Kj and Q on bandwidth and phase margin, which are 
missing in Table 1 due to the order-reduction. 

The parameters of the control loop used in the simulations are listed in Table 2, where V p is the 
biased voltage on the movable mass of the resonator and V mi d is the biased voltage on the drive comb 
fingers. The difference of V p and V mi d determines the Vd c in Equation (4). To compare the simulation 
results with experiment results, the resonance frequency co x and quality Q are especially chosen the 
same as ones in experiments. 



Table 2. Parameters of the control loop simulated in Matlab. 



Parameters (unit) 


Value 


Parameters (unit) 


Value 


K F/V 2 (N/V 2 ) 


3.7e-8 


co x (rad/s) 


8063 x 2ti 


K C /x (F/m) 


1.2e-8 


k (N/m) 


29.6 


K r/c (V/F) 


2.4el2 


Q 


1,368 


Kvga (/V) 


20 


coi p f(rad/s) 


100 x 2tt 


Vref(V) 


0.2 


K P 


5 


V P (V) 


12 


K, 


200 


V mid (V) 


2.5 


X 


1 



The numerical simulations are carried out in Matlab. The simulation results of the transfer function 
in Equation (4) are shown in Figures 4-8a with varying Kp, coi p fi K h and Q. The corresponding step 
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responses at the LPF output terminal of the proposed closed-loop linear model are also given in 
Figures 4-8b. To verify the validity of the proposed linear model, a nonlinear model is built by 
Simulink and step responses at the LPF output terminal of corresponding nonlinear model are shown 
in Figures 4-8c. As a lst-order low pass filter is used in the nonlinear model, the curves in Figures 4-8c 
contain residual resonance-frequency related components and seem noisy. 

Figure 4. (a) Bode diagrams of the control loop (b) step responses of the linear model 
(c) and step responses of nonlinear model varying with the proportional term K P . 
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Figure 5. (a) Bode diagrams of the control loop (b) step responses of the linear model 
(c) and step responses of nonlinear model varying with the cut-off frequency coi p f. 
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Figure 5. Cont. 
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Figure 6. (a) Bode diagrams of the control loop (b) step responses of the linear model 
(c) and step responses of nonlinear model varying with the integral term K h 
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Figure 7. (a) Bode diagrams of the control loop (b) step responses of the linear model 
(c) and step responses of nonlinear model varying with the quality factor Q. 
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Figure 8. (a) Bode diagrams of the control loop (b) step responses of the linear model 
(c) and step responses of nonlinear model varying with the VGA gain coefficient K V ga- 
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Figure 8. Cont. 
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The validity of the order-reduction is proved in Figure 4. In Figure 4a, the variation trend of the 
phase margin first increases and then decreases when K P is doubled. It is noticed that the best phase 
margin is achieved when the zero-pole doublet is created at K P = 10, referring to the resonator 
parameters co x and Q in Table 2. This result indicates that the optimized stability of the control loop is 
achieved when the zero-pole doublet is created, and therefore demonstrates that the validity of the 
order-reduction. The smallest overshoot in Figure 4b at K P = 10 also verifies the conclusion. 

Figure 6a shows that increasing the proportional term Ki is an effective way to enlarge the loop 
gain. However, it causes an obvious decrease in the phase margin and only a little increment in 
bandwidth. The trends of curves in Figure 7a are identical to those in Figure 6a, but the phase margin 
and bandwidth variation quantity is less than those in Figure 6a, which indicates that the stability and 
recovery speed of the control loop is more robust under the variation of quality factor Q than that of 
Kj. It is believed that the phenomenon is decided by zero and poles locations. This conclusion is 
verified by simulations with different parameter settings, but not shown as the long paragraph. The 
curve trends in other figures all agreed with the results in Table 2. 

The step responses of the nonlinear model in Figures 4-8c are similar to those of the linear model in 
Figures 4-8b. However, the simulation responses of the nonlinear model have shorter rise time and larger 
overshoot than those of the linear model. The difference is probably caused by fixing VGA input as its 
final value V re /K rect , which dismisses the action of drive loop in the actual nonlinear model. 

Combining the theoretical and simulation results with the specification optimization order presented 
above, an optimization methodology of system parameters selection can be summarized as follows: to 
achieve a robust, stable and swift control loop, it is better to increase Kj or K VGA to achieve a large loop 
gain and good disturbance rejection capability first, and then pursue an optimized stability by creating 
a zero-pole doublet through adjusting K P . 

3. Circuit Design and Experiments 

To verify the proposed optimization method on a practical gyroscope prototype, AGC-PI based 
closed-loop drive circuits with flexible system parameters adjustment are designed. The block diagram 
is presented in Figure 9. In the drive circuits, an adjustable gain stage is designed to match the 
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mechanical gain variation of different resonators. The PI controller adopts a Gm-C structure, which 
provides convenient and sufficiently tunable properties of K P and Ki [17]. Two high pass filters (HPFs) 
are designed for DC offset cancellation. 

Figure 9. Block diagram of the implemented AGC-PI based closed-loop drive circuits. 
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The drive circuits are implemented in a 0.35 |im CMOS technology. The supply voltage is 5 V. A 
symmetrical doubly decoupled z-axis gyroscope is tested, which adopts an electrostatical drive and 
capacitive sense. The gyroscope is vacuum-packaged and mounted on a printed circuit board (PCB) 
with the implemented drive circuits chip and other necessary external resistances and capacitances, as 
shown in Figure 10. The resonant frequency and the quality factor of the resonator are measured as 
8.063 kHz and 1,368, respectively. 



Figure 10. Test board of the gyroscope prototype. 




To verify the proposed optimization methodology accurately, the velocity-controlled closed-loop 
should be tested working as a linear system during the experiments. Carefully investigating the linear 
model, the input of the VGA module is approximated as the final steady value V re /(2/n). However, in 
the start-up period, the input of the VGA is nearly zero and the gain of VGA is fixed at the maximum 
value due to the overlarge error signal fed into the PI controller. As a result, the velocity-controlled 
closed-loop is broken at the VGA module and the analysis on the linear model above fails. Therefore, 
snatching start-up waveforms as described in [7] cannot be used in our work. 
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In our experiments, a 10 mV step-up signal is applied on the target voltage F^/when the control 
loop is in a steady state. As the step is very small, the VGA input signal can be considered as a 
constant and the velocity-controlled closed-loop is believed to work as a linear system. The LPF's 
output waveform, which reflects the vibration velocity of the resonator, is then recorded with a varied 
proportional term K P and integral term K h K P and Ki in the PI controller are decreased gradually 
by changing the off-chip resistances and capacitances [14]. Figure 11 shows the measured step-up 
waveforms of the LPF's output varied with K P . As K P is increased from 2.5 to 10, the overshoot is 
smoothed gradually, and the rising time and the setting time are improved. When K P is set as 20, the 
overshoot is increased as expected. The test results prove that an optimized stability of the control loop 
is achieved at situation of K P = 10 and Kj = 200, when a zero-pole doublet is almost formed, referring 
to the tested resonators' parameters above. This is exactly the same as the simulation result. The tested 
rising times are given in Table 3 and compared with the simulation results of both linear and nonlinear 
models in Figure 4b,c. It is seen that the simulated nonlinear model's rising times are identical with the 
tested results. As expected, due to the approximation in VGA during the linearization, the rising times 
of linear model show some "delay" compared with the test results. However, their variation trends are 
the same, which also verifies the conclusions in Table 1 in the time domain. Figure 12 shows the 
measured step-up waveforms of the LPF's output varied with Ki and the measurement and simulation 
results of rising time are shown in Table 4. It can be seen that decreasing Kj improves the overshoot of 
the control loop, but increases the rising time and decreases the settling time, which are also in accord 
with the simulation results. 

Figure 11. Measured step-up waveforms of the LPF's output with different proportional 
terms K P . 
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Table 3. Comparison of simulation and measurement rise times with different proportional term K P . 
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Table 4. Comparison of simulation and measurement rise times with different proportional 
term Kj. 
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Disturbance rejection capability (D.R.C.) is an important performance specification of the control 
loop, which is inversely proportional to the loop gain of the negative feedback control loop. To 
evaluate the D.R.C. of the control loop against system parameter variation, experiments are carried out 
by varying the bias voltage V p on the mass. Because the voltage Vdc and the electrical-mechanical 
conversion parameters K V /c in Equation (6) are related to V P9 changing the value of V p can imitate the 
system parameter variation in the control loop. As the TIA's output signal's amplitude, which reflects 
the resonator vibration velocity, is controlled by the negative feedback loop, the variations of Vdc and 
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K v /c will be projected on the amplitude of the drive signal. As a result, the D.R.C. of the control loop 
can be calculated as: 



D.R.C(dB) = 201og 10 



V.P. 



Dr _ amp 



V P 

\ 'TIA_amp J 



(11) 



where V.P.Dr amp and V.P.riAamp are the variation percentage of drive signal's magnitude and TIA' 
output signal's magnitude under different V p values. 

Both the magnitude of the drive signal and the TIA's output are recorded in Table 5 with three 
different values of Kj and under different V p ranging from 16 V to 22 V. The percentages of variation 
of these two signals are calculated by setting V p = 19 V as reference and D.R.C. s are calculated from 
Equation (1 1). As shown in Table 5, the D.R.C. is increased 6 dB when Ki is doubled from 100 to 200. 
This result shows a close agreement with both the theoretical and simulation results. When Kj is 
increased from 200 to 400, the D.R.C. is increased only about 3.6 dB. As the minimum stable bit is the 
sub-mV bit in the AC magnitude measurement, it is believed that the improvement of D.R.C. is limited 
by the noise in the control loop. Figure 13 shows normalized magnitude of drive signal and TIA's 
output with different Kj by varying V p . The benefits of the control loop's D.R.C. can be observed 
directly in Figure 13. 

Table 5. Measured D.R.C of the velocity-controlled closed-loop with different Kj by varying V p . 



Cond. 


V p (V) 


16 


17 


18 


19* 


20 


21 


22 


V.P. (%) 


D.R.C. (dB) 


K T = 100, 


Mag. of Drive 
Sig. (mV) 


419.8 


372.0 


332.8 


299.8 


272.6 


249.6 


229.8 


63.4 


36.2 


K P = 5 


Mag. of TIA's 
Outp. (mV) 


60.4 


60.5 


60.6 


60.7 


60.8 


60.9 


61.0 


0.99 




Ki = 200, 


Mag. of Drive 
Sig. (mV) 


416.0 


368.0 


329.0 


296.8 


269.6 


246.2 


227.0 


63.7 


42.1 


K P = 5 


Mag. of TIA's 
Outp. (mV) 


59.8 


59.9 


59.9 


60.0 


60.0 


60.1 


60.1 


0.5 




Ki = 400, 


Mag. of Drive 
Sig. (mV) 


416.6 


369.0 


330.4 


297.4 


270.4 


247.6 


227.8 


63.5 


45.7 


K P = 5 


Mag. of TIA's 
Outp. (mV) 


59.9 


60.0 


60.0 


60.0 


60.1 


60.1 


60.1 


0.33 





* V.P. Variation Percentage (set V p = 19 V as reference). 

Normalized magnitude of drive signal and TIA's output with different Kj by 



Figure 13. 

varying V p . 




- KI = 100 

- KI = 200 

- KI = 400 




'16 17 18 19 20 21 22 

Bias Voltage Vp/V 



16 17 18 19 20 21 22 
Bias Voltage Vp/V 
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4. Conclusions 

This paper presents a linearization design approach of the 3rd order velocity-controlled closed-loop 
for MEMS vibration gyroscope. As benefits of the built linear model, tedious stability analysis in the 
time domain can be transferred to the frequency domain and the design process of the control loop is 
simplified. The stability criterion and the mathematical expression of each performance specification 
of the 3rd order control loop are presented in this paper. Moreover, an optimization methodology 
which reveals that a robust, stable and swift control loop can be achieved by carefully selecting the 
system parameters following a priority order is summarized. Experiments carried out on a gyroscope 
prototype verify the optimization methodology that an optimized stability of the control loop can be 
achieved by canceling the pole of the resonator with the zero in the PI controller, and disturbance 
rejection capability (D.R.C) of the control loop can be improved by increasing the integral term Kj. It 
is also found that noise in the control loop limits the D.R.C. of the control loop. As a result, the noise 
performance of the control loop will be studied in the future work. 

Acknowledgments 

The authors gratefully acknowledge National Natural Science Foundation of China (project 
61106025) and the CAS/SAFEA International Partnership Program for Creative Research Teams for 
financial support. 

Conflicts of Interest 

The authors declare no conflicts of interest. 
References 

1. Yazdi, N.; Ayazi, F.; Najafi, K. Micromachined Inertial Sensors. Proc. IEEE 1998, 86, 1640-1659. 

2. Nasiri, S. A critical review of MEMS gyroscopes technology and commercialization status. 
Avaliable online: http://invensense.com/mems/gyro/documents/whitepapers/MEMSGyroComp.pdf 
(accessed on 25 May 2012). 

3. Nasiri, S. A critical review of the market status and industry challenges of producing consumer 
grade MEMS gyroscopes. Avaliable online: http://invensense.com/mems/gyro/documents/whitepapers 
(accessed on 25 May 2012). 

4. Shaeffer, D.K. MEMS inertial sensors: A tutorial overview. IEEE Commun. Mag. 2013, 57, 100-109. 

5. Sung, W.-T.; Sung, S.; Lee, J.-Y.; Kang, T.; Lee, Y.J.; Lee, J.G. Development of a lateral 
velocity-controlled MEMS vibratory gyroscope and its performance test. J. Micromech. 
Microeng. 2008, 18, 055028, doi:10.1088/0960-1317/18/5/055028. 

6. Xia, D.; Chen, S.; Wang, S. Development of a prototype miniature silicon microgyroscope. 
Sensors 2009, 9, 4586^605. 

7. Cui, J.; Chi, X.Z.; Ding, H.T.; Lin, L.T.; Yang, Z.C.; Yan, G.Z. Transient response and stability of 
the AGC-PI closed-loop controlled MEMS vibratory gyroscopes. J. Micromech. Microeng. 2009, 
19, 125015, doi:10.1088/0960-1317/19/12/125015. 



Sensors 2013, 13 



12580 



8. Aaltonen, L.; Halonen, K.A.I. An analog driving loop for a capacitive MEMS gyroscope. Analog 
Integr. Circuits Signal Process. 2010, 63, 465-476. 

9. Chen, Y.-C; Closkey, R.T.M. A control and signal processing integrated circuit for the 
JPL-Boeing micromachined gyroscopes. IEEE Trans. Control Syst. Technol. 2005, 13, 286-300. 

10. Xia, D.; Yu, C; Wang, Y. A digitalized silicon microgyroscope based on Embedded FPGA. 
Sensors 2012, 12, 13150-13166. 

11. Aaltonen, L.; Kalanti, A.; Pulkkinen, M; Paavola, M; Halonen, K.A.I. A 2.2mA 4.3 mm 2 ASIC 
for a 1000 deg/s 2-Axis Capacitive Micro-Gyroscope. IEEE J. Solid-State Circuit 2011, 46, 
1682-1692. 

12. Cui, J.; Giio, Z.Y.; Yang, Z.C.; Hao, Y.L; Yan, G.Z. Electrical coupling suppression and transient 
response improvement for a microgyroscope using ascending frequency drive with a 2-DOF PID 
controller. J. Micromech. Microeng. 2011, 21, 095020, doi: 10. 1088/0960-13 17/21/9/095020. 

13. Closkey, R.T.M. ; Vakakis, A. Analysis of a Microsensor Automatic Gain Control Loop. In 
Proceedings of American Control Conference, San Diego, CA, USA, 2 June 1999; 
pp. 3307-3311. 

14. Loveday, P.W.; Rogers, CA. The influence of control system design on the performance of 
vibratory gyroscopes. J. Sound Vib. 2002, 255, 417-432. 

15. Leland, R.P. Adaptive control of a MEMS gyroscope using lyapunov methods. IEEE Trans. 
Control Syst. Technol. 2006, 14, 278-283. 

16. Xia, D.; Chen, S.; Wang, S.; Li, H. Microgyroscope temperature effects and compensation-control 
methods. Sensors 2009, 9, 8349-8376. 

17. Wu, H.; Yang, H.; Cheng, X.; Yin, T. Integrated Gm-C based PI Controller for MEMS Gyroscope 
Driveloop. In Proceddings of IEEE 9th International Conference on ASIC, Xiamen, China, 
25-28 October 201 1; pp. 925-928. 



© 2013 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article 
distributed under the terms and conditions of the Creative Commons Attribution license 
(http ://creati vecommons . org/licenses/by/ 3.0/). 



